Double-periodic quasi-periodic graphene superlattice: 
non-Bragg band gap and electronic transport 



1 Department of Physics, Shanghai University, 200444 Shanghai, China 

2 Departamento de Qmmica-Ffsica, UPV-EHU, Apdo 644, 48080 Bilbao, Spain 

Abstract. Electronic band gap and transport in quasi-periodic graphene superlattice of 
double-periodic sequence have been investigated. It is found that such quasi-periodic structure 
can possess a zero-averaged wave number (zem-k) gap which associated with an unusual Dirac 
point. Different from Bragg gap, the zeio-k gap is less sensitive to the incidence angle, and 
robust against the lattice constants. The locations of Dirac point and multi-Dirac -points in 
the graphene superlattices of various sequences are also compared. The control of electron 
transport over the zero-fc band gap in graphene superlattice may facilitate the development of 
many graphene-based electronics. 




Submitted to: J. Phys. D: Appl. Phys. 



$ Corresponding author, xchen@shu.edu.cn 



Double-periodic quasi-periodic graphene superlattice.. 



2 



1. Introduction 

Since the discovery of graphene in 2004, electronic band gap and transport in graphene 
have attracted considerable attention because of the intriguing physics as well applications 
in graphene-based nanoelectronics E El El SJ. Graphene has a unique band structure with 
the conductance and valance bands touching at Dirac point, which leads to many unusual 
properties of transport [4], like half-integer quantum Hall effect, Klein tunneling, and the 
minimal conductance. So far, the transport properties including Klein tunneling and resonant 
tunneling have been extensively investigated in various graphene-based heterostructures 
ranging from single barrier jH [6l to superlattice with electrostatic potential and magnetic 

barriers £Zl El |5l HD H21 H31 US H31 - 

Several works have been devoted to the new Dirac point in the band structure of graphene 
superlattice (GSL) 111 [12 [HE EL smce sucn new DP, associated with zero-averaged wave 
number (zero-fc) gap, is of benefit to the controllability of electronic transport. Different 
from Bragg gaps, the zero-fc gap emerged with new Dirac point is robust against the lattice 
constant, structural disorder [fT2l . and external magnetic field [13J. Besides, the structure 
ordering of quasi-periodic GSL is between periodic and disordered systems, which could give 
rise to many interesting and significant phenomena, for example, fractal spectrum and self- 
similar behavior lfT8l[T9l . Recently, we have studied the zero-fc gap and electronic transport in 
Fibonacci quasi-periodic GSL |[T6ll . Shortly afterwards, Ma et al. IfTTll found that the zero-fc 
gap also exists in other quasi-periodic GSL of Thus-Morse sequence. 

As we know, the main examples of quasi-periodic models are the Fibonacci (FB), Thue- 
Morse (TM) and Double-Periodic (DP) structures [|20l . These quasi-periodic superlattice can 
be generated by the following substitution rules: A — > AB and B — > A for FB, A — * AB and 
B — > BA for TM, and A — » AB and B — > AA for DP, where A and B are two different barriers 
in the quasi-periodic superlattice. The DP quasi-periodic superlattice is obviously different 
from FB and TM ones. But the number of barrier layers in DP sequence increases as N = 2" 
(like the TM sequence), where n indicates the iteration order, and barriers B appear always 
isolated (like the FB sequence) [20]. The electronic properties of FB, TM and DP lattice are 
of significance [f2TTl . So it is worthwhile to study the electronic band gap and transport in DP 
GSL and make comparisons with different quasi-periodic systems. Combing the results in FB 
and TM GSL lfT6l[P71 . the controllable conductance and shot noise in DP GSL may result in 
the practical applications of electron transport in quasi-periodic GSL. 

2. Theoretical Model 

Consider a graphene-based DP quasi-periodic superlattice which is arranged by DP 
substitution rule, A — » AB and B — » AA, where A and B are two different barriers, the 
parameters are denoted by the potential barrier height V A and width d A , the potential well 
height Vb and width J B . For an n-fh DP sequence, S„, it contains elements A and B, and 
follows the relation, S„ = S n -iS n _ v Si = S n -iS n -i (for n > 1) with Si = AB and S\ = AA, 
which leads to S 3 = ABAA, S 3 = ABAAABAB, S 4 = ABAAABABABAAABAA and so on. 
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Figure 1. (Color online) Trace-maps for DP quasi-periodic GSL, (S„)' n , of variable order n 
under (a) d A — ds — 20 nm, (b) cIa - 2ds - 30 nm, where other parameters are Va = 50 meV, 
Vb — meV, 6q - 20°, and m - 12. The horizontal dashed red line denotes the new Dirac 
point's location. 



The Hamiltonian of change carriers near the K point inside a monolayer graphene is given 
by, H = -ifivpcr • V + V(x), where the Fermi velocity v F ~ 10 m/s, and cr = (a x , a } ) are the 
Pauli matrices and V(x) denotes the potential barrier or well. The solution of H, acting on the 
electronic pseudospin wave functions, results in the transfer matrix [|T2l : 



' cos,(qjAx-8j) .s'm(qjAx) 



■km / a p i \ cos 6j cos 6j 

M j\L\X, tL , Ky) — .smiqjhx) co$(qjAx+9j) 

y ' COS0j cosBj 



(1) 



which connects the wave functions at x and x + Ax inside the jth potential with 8j = 
axcs'm(ky/kj), where kj = (E - Vj)/~hv F , k y and qj are the y and x components of wavevector, 
qj = sign(k j)(k 2 j - k 2 ) 112 for kj > k 2 „ otherwise qj = i(k 2 - k 2 ) i/2 . Consequently, the 
transmission coefficient t = t(E, k y ) is found to be 

t = ^2 , (2) 

(m 2 2e~' e ° + m n e l0 <) - m l2 e' {6 >~ eo) - m 2 i ' 

where 6 and 6 t are incidence and exit angles, and m !; (z, j = 1,2) is the matrix element 
of total transfer matrix, X[S n ] = \\ N j =l Mj{dj,E,k y ). Once the transmission coefficient is 
obtained, we can calculate the total conductance G [221 and the Fano factor F 11231 in terms 



of G = Go £ 12 T cos OodOo, and F = JT^ T(l - T) cos 6 d9 / ^ T cos 6 d9 , where T = \t\ 2 
and Go = 2e 2 mv F L y /h 2 , with L y is the width of the graphene stripe in the y direction. 

In Fig. ED (a) and (b), we plot the trace maps for quasi-periodic GSL of DP sequences, 
S™, where (a) d A = d B = 20 nm, (b) d A = 2d B = 30 nm, the incidence angle 6 = 20°, V A = 50 
meV, Vb = meV, and m = 12. Similar to the results for FB and TM sequences lfT6l [TTll . 
we find that there are several broad forbidden gaps open for each DP level, and the passing 
bands are split into more and more subbands as DP sequence order n increases. Among these 
forbidden gaps, the center position of the zero-fc gaps, which are denoted by a dashed red line, 
are almost the same for different DP sequences in Fig. ED (a) and (b). These zero-fc gaps are 
different from Bragg gaps, and thus are robust against the lattice parameters. 
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Figure 2. (Color online) Band structures for DP quasi-periodic GSL, (ABAA) 12 , where (a) 
oU = ds = 20 nm, (b) <Ia - 2ds - 30 nm, and the other parameters are the same as those in 
Fig. Q] The horizontal dashed red line denotes the new Dirac point's location. 




Figure 3. (Color online) Transmission spectra for DP quasi-periodic GSL, (ABAA) , where 
(a): cLa/cLe — 1, flU = 15 nm (dashed black line), <Ia - 20 nm (solid red line), (b) cLa/cLe = 2, 
d-A = 20 nm (dashed black line), oU = 30 nm (solid red line), and the other parameters are the 
same as those in Fig. [TJ 



To understand better, we further calculate the electronic dispersion, for example the 
second DP sequence (S 2 = ABAA), at any incidence angle 

cos(£,A 2 ) = l -Tr{M A M B M A M A ), (3) 

where B x is the x component of Bloch wave vector, and A 2 = 3d A + dg is the length of the 
unit cell. From Eq. ©, we can finally get the electronic band structures, as shown in Fig. [2] 
(a) and (b), which correspond to (a) d A = d B = 20 nm, (b) d A = 2d B = 30 nm. In general, the 
location of new Dirac point inside the zero-k gap is determined by 

N N 

k = Y J kjdj/Yjd J = 0, (4) 



Double-periodic quasi-periodic graphene superlattice.. 



5 




E (meV) 



Figure 4. (Color online) Transmission spectra for DP quasi-periodic GSL, (ABAA) i , where 
(a): d-A — ds — 20 nm, (b): d A = 30 nm, ^ = 15 nm, 6q = 15° (dashed black line), 6q = 20° 
(solid red line), 9q = 25° (dotted magenta line), and the other parameters are the same as those 
in Fig. Q] 



which is valid for both periodic and aperiodic GSLs [T2l [T6l [TVJ . Thus we achieve the new 
Dirac point located at 

V A + V B ■ d B l(T n d A ) , c , 

E n = , (5) 

1 + d B /(T n d A ) 

which depends on the ratio of numbers of layer A and B, that is, r„ = N A /N B . Here we would 
like to emphasize that the location of Dirac point is generally dependent of the order n. But 
when TM sequence is considered, the location of Dirac point does not depend on the order n 
because of t„ = 1 IfTTTl . The different behaviors will be compared later in the different GSLs of 
FB, TM and DP sequences. Besides, we see the Dirac point also depends on the ratio d A /d B . 
As shown in From Fig. [2](a) and (b), the position of Dirac point is shifted from E = 37.5 meV 
for d A /d B = 1 to E = 42.875 meV for d A /d B = 2. Furthermore, Fig. [3]shows the transmission 
spectra for DP quasi-periodic GSL, (ABAA) 12 . It is demonstrated that the zero-k gap is robust 
against the ratio of lattice parameters, d A /d B , while other Bragg gaps are sensitive to lattice 
parameters. From the point of view of application, this provides the flexibility to control the 
electron transport by adjusting the ratio of lattice parameters, d A /d B , in the various sequences. 

Next, the influences of the incidence angle, 8 Q , on the transmission spectra are also 
studied. Obviously, the zero-k gap associated with Dirac point is insensitive to the incidence 
angle, as shown in Fig. HI On the contrary, the Bragg gap shifts upward in energy as the 
incidence angle increases. Compared to the weak dependence of zew-k gap on the incidence 
angles, the Bragg gaps change sensitively with respect to incidence angle 6>o- 

Additionally, the multi-Dirac -points could appear in the GSL with periodic |T2l and 
aperiodic [fTTl structures. Fig. ^demonstrates that the extra Dirac points located at k y ± 0, 
can also emerge in the GSL of DP sequence. We show the energy band gaps versus the 
lattice constants in the case of d A = d B = d. The zero-fc gap, associated with the multi- 
Dirac -points, oscillates with changing the lattice constants, thus open and close periodically, 
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Figure 5. (Color online) Electronic band gaps depending on the lattice constants d& — ds — d 
with k y = 0.01 nm -1 , where (a) S 2 = ABAA for DP sequence, (b) S 2 = ABBA for TM 
sequence, and the other parameters are the same as those in Fig. [2] The horizontal dashed red 
line denotes the Dirac point's location for k y = 0, where (a) E = 37.5 meV and (b) E — 25 
meV. 

while the other Bragg gaps are significantly shifted. However, the multi-Dirac-points in DP 
sequence are different from those in TM sequence. To compare them, we choose S % = ABAA 
for DP and S 2 = ABBA for TM sequence. We see from Fig. [5] (a) that the multi-Dirac-points 
appear at the same energy, E = 36.4 meV, which is different from the location, E = 37.5 meV, 
of Dirac point for k y = 0. Moreover, the energy corresponding the multi-Dirac-points in DP 
sequence will be close to E = 37.5 meV, when k y — » 0. However, the Dirac point for k y = 
and multi-Dirac-points for k y ± happen at the same energy E = 25 meV in TM sequence 
with t„ = 1 . It tells us that though the total numbers of layer A and B are the same for TM and 
DP sequences, the behaviors of the multi-Dirac-points in TM and DP sequences are totally 
different, because the ratio, t„ = N A /N B , are not equal. 

Regarding the electronic transport, the total conductance and the Fano factor are also 
interesting. In Fig. [51 we calculate the total conductance G/Gq and the Fano factor F in 
quasi-periodic GSL of DP squence, (ABAA) 12 . Again, we have found serval unique features 
associated with new Dirac point and zero-£ gap, as shown in Fig. [6] Firstly, the curve of 
angular-averaged conductance reaches its minimum at the Dirac point, and forms a liner cone 
around Dirac point. Meanwhile the Fano factor reaches the value of 1/3 approximately [|23l . 
Secondly, the conductance and the Fano factor will shift with the ratio of Ja/^b- Actually, 
this will be useful to modulate the conductance of such GSL by changing the ratio of lattice 
constants. 

Finally, what we should mention is that the transmission spectrum and corresponding 
electronic transport properties including conductance and the Fano factor in the DP sequence 
are quite different from those in FB and TM sequences. To clarify it, we shall compare the 
Dirac point in three different sequences. In Fig. [71 we present the location of Dirac point 
versus the variable order n in FB, TM and DP sequences. In TM sequences, the location of 
new Dirac point does not depend on the order n, because, the numbers of layer A and B are 
the same as mentioned above. As for FB and DP, the numbers of layer A and B are different 
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Figure 6. (Color online) Conductance G/Go [(a) and (b)] and Fano factor F [(c) and (d)] 
versus Fermi energy E in DP quasi-periodic GSL, (ABAA) n , where (a,c) c/a/^b = 1, 4a = 15 
nm, (b, d) d A /d B = 2,cLa = 20 nm, and the other parameters are the same as those in Fig. \T\ 




Figure 7. (Color online) Location of Dirac point versus variable order n in FB, TM and DP 
sequences, where (a) dA/dB = 1, (b) dA/ds - 2, and the other parameters are the same as those 
in Fig. Q] 



in each order n (n > 1). As a consequence, the location of Dirac point depends on the variable 
order n, but it will becomes stable for n » 1, since the ratio of numbers of layer A and B, 
N A /N B , trends to t„ = (1 + V5)/2 for FB sequence OH HI and T n = 2 roughly for DP 
sequence ||20l . These difference presented here are quite natural, because the physical origin 
of zero-fc gap and new Dirac point is the total zero phase, suggested by Eq. ©. 



3. Summary 

To summarize, we have studied the electronic band gap and transport in GSL of DP sequence 
using the transfer matrix method. We have found that the Dirac point, multi-Dirac -points and 
associated zero-k gaps have emerged in such quasi-periodic structure. The zero-fc gap and 
Dirac point are robust against the lattice constants and incidence angle. The results are also 
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compared with those in GSL of FB and TM sequences. We hope that all these results may 
lead to applications in the control of electron transport in GSL. 
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